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基于 变异 交叉 方程 与 进化 选择 机 制 的 回溯 优化 改进 算法 * 
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摘 要 : 针对 回溯 搜索 优化 算法 存在 的 收敛 速度 慢 ， 容 易 陷入 局 部 最 优等 问题 ， 提 出 了 一 种 改进 算法 。 首 先 利 用 t 分 
布 产 生变 蜡 尺 度 系 数 ， 加 快 了 算法 收敛 速度 ; 接着 完善 交叉 方程 结构 ， 引 入 最 优 个 体 控制 种 群 搜索 方向 ， 有 效 提高 了 
算法 开发 能 力 ; 最 后 提出 进化 选择 机 制 ， 引 入 差分 进化 算法 变异 因子 ， 一 定 概率 下 以 较 差 解 替 换 较 优 解 ， 训 免 算 法 陷 
入 局 部 最 优 。 在 数值 实验 中 ， 选 取 了 15 个 测试 泡 数 进行 仿真 测试 ， 并 与 5 种 表现 良好 的 算法 进行 了 比较 ， 结 果 表明 ， 
该 算法 在 收敛 速度 及 搜索 精度 方面 有 明显 优势 。 
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Backtracking search optimization algorithm based on mutation and 
crossing equations and evolutionary selection mechanism 


Zhao Linjing!'™ 't, Ge Baozhenla 了 ,Chen Leila 1b, ?1 
(1. a. School of Precision Instruments & Opto-Electronics Engineering, b. Key Laboratory of Opto-Electronics Information & 
Technical Science of Ministry of Education, Tianjin University, Tianjin 300072, China; 2. School of Information Engineering, 
Tianjin University of Commerce, Tianjin 300134, China) 


Abstract: According to the slow convergence and easiness to trap in local optimum of backtracking search optimization 
algorithm, this paper presented an improved algorithm. The method proposed a mutation scale factor based on t distribution 
firstly to speed up the convergence rate; then the algorithm improved the structure of crossover equation and introduced the 
optimal individual to control the direction of population search, which effectively improved the development capability. Finally, 
the algorithm proposed the evolutionary selection mechanism, introduced the mutation factor of differential evolution algorithm 
and replaced the optimal solution with worse solution under a certain probability, which can avoid algorithm to fall into the local 
optimum. The numerical experiments selected 15 test functions for simulation and compared with 5 well-behaved algorithms. 
The results show that the proposed algorithm has obvious advantages in terms of convergence rate and search accuracy. 
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引言 人 吕 提 出 的 蜂 群 改进 算法 GABC， 李 文 斌 等 人 中 提出 的 差分 进 
呈 化 改进 算法 CNS-DE， 姜 凤 利 等 人 叫 提 出 的 粒子 群 优化 改进 算 
近 些 年 ， 进 化 算法 以 其 结构 简单 、 参 数 较 少 、 易 于 实现 等 法 IPSO， 以 及 Meng 等 人 四 提出 的 鸡 群 算法 (chicken swarm 


寺 点 受到 广泛 应 用 ， 各 种 进化 算法 层 现 铸 出 ， 如 经 典 的 遗传 算 。 ”optimization, CSO) ，Zheng 等 人 中 提出 的 水 波 算法 (water wave 
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法 (genetic algorithm, GA) 中 、 粒 子 群 优化 算法 (particle swarm optimization, WWO) ， 高 江 锦 等 人 09 提 出 的 免疫 否定 进化 选择 
ptimization, PSO ) 外、 差分 进化 算法 (deferential evolution, DE) 算法 (immune evolution negative selection algorithm,IENSA)〉 等 
、 人 工蜂 群 算 法 (artificial bee colony algorithm, ABC)》 也 等。 新 型 进化 算法 。 
然而 随 着 高 维 多 模 态 复杂 优化 问题 出 现 的 愈加 频繁 ， 使 全 局 优 回溯 搜索 优化 算法 (backtracking search algorithm, BSA ) 
化 算法 面临 着 更 大 的 挑战 。 对 此 ， 学 者 们 在 现 有 进化 算法 的 基 Civicioglu00 于 2013 年 提出 的 一 种 新 型 进化 算法 。 算 法 框架 与 
础 上 ， 提 出 了 大 量 改进 算法 以 及 新 的 优化 算法 。 例 如 ，Zhu 等 差分 进化 算法 类 似 ， 但 县 体操 作 有 本 质 区 别 053。BSA 采用 
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种 群 控制 搜索 方向 ， 加 之 其 高 效 新 颖 的 混合 交叉 策略 ， 使 其 各 
项 性 能 都 取得 了 不 错 的 成 绩 。 


Chin 正四 天 ， 
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的 随机 数 ; mixrate 为 交叉 概率 ; mixrate :rand:D 为 交叉 长 度 ; 
& = randperm(D) ， 为 [1,2:…D] 重新 排序 后 的 整数 向 量 。 BSA 通过 


但 整体 来 看 ，BSA 探索 能 力 较 强 而 开发 能 力 有 待 提 高 ， 导 


致 收敛 速度 较 慢 。 针 对 这 一 缺点 ， 本 文 提 出 了 基于 t 分 布 的 变 
异 系数 FE， 以 及 加 入 最 优 个 体 引导 的 新 交叉 方程 ， 在 保证 探索 
能 力 的 基础 上 ， 大 大 提高 了 算法 的 收敛 速度 。 同 时 ， 提 出 进化 
选择 机 制 , 引入 差分 进化 算法 变异 算 子 DE/rand/1, 增加 全 局 探 
索 能 力 ， 避 和 免 算法 早熟 收敛 。 最 后 ， 选 取 15 个 测试 函数 ， 对 包 
括 本 文 改进 算法 在 内 的 六 种 算法 进行 数值 实验 ， 结 果 表 明 ， 本 
文 改进 算法 显著 提高 了 BSA 的 收敛 速度 和 精度 。 


1 ”回溯 搜索 优化 算法 介绍 


BSA 是 一 种 基于 群体 智能 的 进化 算法 , 算法 结构 与 其 他 进 
化 算法 类 似 , 主要 包括 五 个 步骤, 分 别 为 种 群 初始 化 、 选 择 I 、 
种 群 变异 、 种 群 交 义 、 选 择 II。 
1) 种 群 初始 化 
BSA 初始 化 种 群 的 方法 为 随机 产生 ,方程 如 式 (1)(2) 所 示 。 
Pop ; ~ U(Low ,Up,) (1) 


Oldpop; ; ~ U (Low,, UP;) (2) 
其 中 ，Pop 为 初始 种 群 ，Oldpop 为 历史 种 群 ，j (1,2...N}， 
je{1,2…D}; 和 是 种 群 个 数 ;，D 是 种 群 维 数 ，Low 和 Up 分 别 


调用 这 两 种 等 概率 的 Map 产生 方式 ， 有 效 控制 新 种 群 7 中 变 
异 元 素 的 个 数 。 
在 新 种 群 工 确定 后 ， 对 种 群 中 的 元 素 进 行 边界 检测 ， 
出 搜索 范围 ， 则 按 式 (1) 产 生 新 的 个 体 。 
5) 选择 开 
通过 贪 禁 法 则 在 新 种 群 与 初始 种 群 中 选择 适应 度 值 较 好 的 
个 体 ， 并 记录 当前 全 局 最 优 解 和 对 应 的 解 向 量 ， 同 时 更 新 初始 
种 群 。 更 新 方程 如 式 (8) 所 示 。 


知 超 


TH 


到 fimess(T:;) < fimess(Pop;) 


Pop; = (8) 


Pop; otherwise 


复 上 述 过 程 ， 直 到 满足 终止 条 件 ， 


[wll 
Th 


最 后 输出 最 优 解 。 


2 “BSA 改进 算法 


BSA 算法 全 局 探索 能 力 较 强 而 局 部 开发 能 力 稍 显 弱势 ， 究 
其 原因 是 该 算法 仅 在 历史 种 群 的 引导 下 进行 变异 ， 缺 少 最 优 个 
体 的 引导 。 针 对 这 一 特点 ,本 文 对 变异 及 交叉 方程 进行 了 改进 ; 


另 一 方面 , 本 文 提 出 了 算法 进化 率 P, 引入 差分 进化 变异 因子 ， 
使 算法 后 期 种 群 多样 化 ， 避 免 算 法 陷入 局 部 最 优 。 本 文 将 所 提 


出 的 改进 算法 简 记 为 MEBSA。 
2.1 变异 尺度 系数 改进 
原始 BSA 算法 的 变异 尺度 系数 严 =3*randn ， 分 布 范围 过 


为 搜索 空间 的 下 界 和 上 界 ; U 为 随机 均匀 分 布 函数 。 
2) 选择 I 
BSA 通过 历史 种 群 控制 搜索 方向 ,选择 I 的 目的 是 在 每 次 


迭代 开始 时 选择 新 的 历史 种 群 。 选 择 策略 如 式 (3) 所 示 。 
if a <b then Oldpop = Pop (3) 


其 中 : a、b 是 (0,1) 上 服从 均匀 分 布 的 随机 数 。 在 确定 OUdpop 


之 后 ， 将 Oldpop 中 种 群 的 位 置 进行 重新 排列 ， 公 式 如 下 : 
Oldpop = randperm(Oldpop) (4) 
其 中 : randperm 是 重新 排序 函数 。 
3) 变异 
BSA 的 变异 方程 如 (5) 所 示 。 
M = Pop + F : (Oldpop -Pop) (5) 


其 中 : FF 为 变异 尺度 系数 用 以 控制 搜索 方向 矩阵 (Oldpop — Pop) 
的 幅度 ，F =3*randn , randn 是 一 服从 标准 正 态 分 布 的 随机 数 。 
4) 交叉 


BSA 提出 了 一 种 新 型 的 交叉 策略 , 通过 设置 混合 比例 参数 
来 控制 种 群 粒子 交叉 个 数 。 有 具体 公式 如 式 (6) 所 示 。 
M.., Map, ;=1 
J ,J 
Dj | J Map; ; =0 (6) 


其 中 :Map ; 为 NXD 的 二 元 整数 矩阵 ;初始 赋值 为 1。 具体 公 
式 如 式 (7) 如 示 。 


=0, a<b 


J ,(ul l:mixrate-rand:D|) (7) 


Map a>b 


i,randi(D) 一 


其 中 : randi(D) 为 从 [0, DP] 中 随机 取 一 个 整数 , rand; a、5 为 [0,1] 


广 ， 容 易 产 生 较 大 或 较 小 的 极 值 ， 不 利于 种 群 的 快速 收 化 。 为 
此 ， 本 文 提 出 了 一 种 基于 t 分 布 的 变异 尺度 系数 。 其 表达 式 为 
F =2* Vtrnd(df) (9) 


其 中 : trnd(q) 是 一 服从 t 分布 的 随机 数 。 
t 分 布 只 有 一 个 参数 , 即 自由 度 df, 其 密度 函数 较 于 标准 正 
态 分 布 的 密度 函数 表现 为 厚 尾 , 且 其 形状 随 着 df 的 增 大 而 无 限 


接近 于 正 态 分 布 。 因 此 ， 相 比 于 正 态 分 布 ， 通 过 自由 度 ij, 五 
以 很 方便 地 控制 t 分 布 的 形状 ， 且 因 其 厚 尾 的 特点 ， 在 整体 分 
布 趋势 较为 集中 的 基础 上 ， 也 可 以 一 定 概率 取 到 处 于 边界 的 极 
值 ， 从 而 避免 算法 过 分 收敛 , 陷入 局 部 最 优 。 本 文 取 d 广 10。 基 
于 分布 的 变异 尺度 系数 与 原始 系数 绝对 值 比较 ， 在 不 同 区 间 
的 简要 分 布 情况 如 表 1 所 示 。 


表 1 各 区 间 简 要 分 布 概率 


区 间 (0,0.5) (0.5,D) (1,1.5) (152) (24) >4 
2* tmnd(df) 0.0487 0.1412 0.2229 0.2447 0.3404 0.0021 
3*randn 0.132 0.128 0.121 0.112 0.323 0.1824 


表 1 可 以 看 出 ， 与 原始 系数 相 比 ， 新 的 变异 系数 产生 的 
扰动 较为 集中 ， 且 不 会 产生 过 大 的 数值 ， 加 快 了 种 群 的 收敛 速 
度 ， 一 定 程度 上 提高 了 算法 的 局 部 开发 能 力 。 

2.2 ”交叉 方程 改进 
通过 对 BSA 算法 的 分 析 可 知 , 该 算法 因 其 独特 的 变异 交叉 
操作 ， 全 局 探索 能 力 强 ， 但 因 其 仅 在 历史 种 群 的 引导 下 控制 搜 
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索 方向 ,导致 该 算法 的 收敛 速度 较 慢 ， 
针对 这 一 缺点 ， 本 文 算法 引入 最 优 个 


的 交叉 方程 ， 公 式 如 式 (10)(11) 所 示 。 


其 中 : 
和 矩阵， 


来 ， 


人 j +C(Gpop; ; 一 Popi 让 

7 LPom +C(Gpopi ;— Pop; ;), Map; ;=0 
Gpop = rempat((normrnd(],1).* Gminimizer), N) 

正 态 分 布 的 随机 数 ;，Gpop 是 一 

当前 种 群 最 优 个体 的 ] 

完整 ， 不仅 


C 为 一 服从 标准 


矩阵 元 素 
BSA 算法 的 交叉 变异 方程 趋 于 


本 的 引导 ， 


Map; j = 


局 部 开发 能 力 还 不 够 强 。 
提出 了 一 个 新 


(10) 


(11) 
个 NxD 


还 有 当前 最 


饮 个 体 控制 搜索 方向 。 


E 态 分 布 组 成 。 这 样 一 
历史 种 群 
值得 一 提 的 是 ， 本 文 没 


的 引 


有 直接 采 


动 ， 
力 。 
2.3 


着 循环 次 数 的 增多 ， 仍 然 面 
的 多 峰 多 极 值 函数 ,存在 可 能 无 法 搜索 到 最 优 人 


用 当前 全 局 最 优 个 体 引导 
使 算法 在 大 大 提高 


进化 选择 机 制 


， 而 是 加 上 了 
了 收敛 速度 之 外 ， 


原始 BSA 算法 在 初期 迭代 中 


本 文 提 出 了 进化 选择 机 制 


代 ， 


个 体 ， 引 入 差分 进化 
个 体 , 并 替换 掉 当 前 
率 将 较 差 的 解 代 蔡 较 优 解 ， 
避免 早熟 收敛 。 新 个 体 产 生 公 式 如 


其 中 : 
随机 数 ; rl#¥r2e {1,2,…N} 
2 =0.5, 


计算 一 次 了 P。 当 


使 


法 中 后 大 
式 (12) 所 示 。 


的 种 群 


Newpop; = Popy(i) +W(Pop,] -Pop,,) 


ie {1,2,…0Q}， 


v(i)=cr(i), cr = randperm(N); 


E 态 分 布 的 扰 
也 保证 了 全 局 探索 能 


有 较 好 的 探索 能 力 , 然而 随 
鱼 着 早熟 收敛 的 问题 ， 对 于 较 复杂 


直 的 缺点 。 因 此 ， 
， 引 入 算法 进化 率 P， 设 定 每 隔 per 
尸 小 于 所 设 阔 值 。 时 ， 随 机 选 
法 变异 算 子 DE/rand/1， 产 生 0 个 新 的 
1 种 群 及 历史 种 群 中 的 相应 个 体 ,以 一 定 的 概 
更 加 多 样 化 ， 


取 个 当前 


(12) 


玉 取 [-1,1] 的 


QO =ceil(0.3*N)。 其 中 : 


出 改进 BSA 算法 的 流程 : 
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~ a 
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16 


17 


种 群 初始 化 
while 不 满足 算法 停止 条 件 do 
选择 I， 选 取 0Ldpop 


进行 
if rand<rand do 


Map =0 


i,(uf :mixrate-rand:D |) 
else 
Map; ,aicpy =0 


end if 


按 式 (16) 进 行 变 异 操作 
按 式 (6) 进 行 交 叉 操作 


if P<6.5 


for i=1: ceil(0.3*N) 


按 式 (12) 产 生 新 的 个 体 ， 并 


end for 


蔡 换 


end if 


进行 选择 工 ， 更 新 Pop 


end while 


ceil 为 向 上 取 整 函数 。 下 面 给 


互 不 相同 。 本文 算法 设 定 per=100， 


给 
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3 ”数值 实验 及 结果 分 析 


为 了 测试 本 文 改进 算法 的 效果 ， 选 取 了 15 个 测试 函数 进 
行 数值 实验 ， 同 时 与 标准 BSA、BSA 改进 算法 COBSAU3]、DS 
09、MABC 算法 05 以 及 CSO 算法 进行 了 对 比 实验 。 表 2 给 出 
了 15 个 测试 函数 的 名 称 、 维 数 、 搜 索 空 间 和 最 优 
秆 0 。 其 中 ，F1~F7 是 简单 的 单 峰 测试 函数 ; F8~F9 是 
多 峰 少 极 值 函 数 ，F10~F15 是 多 峰 多 极 值 函数 。 这 些 函数 的 局 
部 极 值 个 数 随 着 函数 维 数 的 增加 呈 指 数 级 增长 ， 这 类 函数 对 大 


多 数 算法 都 是 很 难 求解 的 。 
表 2 15 个 函数 的 简单 描述 
函数 名 称 维 数 搜索 空间 最 优 值 
Fl Sphere 60 [-100,100] 0 
F2 Schwefel 2.22 60 [-10,10] 0 
F3 Sumsquare 60 [-10,10] 0 
F4 Exponential 60 [-1.28,1.28] 0 
FS Schwefel 1.2 10 [-30,30] 0 
F6 Elliptic 30 [-100,100] 0 
F7 R-H-Ellipsoid 80 [-100,100] 0 
F8 Zakharov 60 [-5,10] 0 
F9 Salomon 60 [-100,100] 0 
F10 Alpine 60 [-10,10] 0 
F11 Rastrigin 100 [-100,100] 0 
F12 Griewanlk 60 [-600,600] 0 
F13 Schaffer 10 [-100,100] 0 
F14 Weierstrass 60 [-0.5,0.5] 0 
F15 NCRastrigin 50 [-5.12,5.12] 0 
3.1 停止 条 件 制定 及 参数 设 定 


Ca) 所 求 结果 与 测试 函数 的 理论 最 优 解 绝 对 值 之 差 小 于 
Matlab 系统 最 小 正 浮 点 数 4.9407e-324， 即 eps (0) 时 ， 停 目 
(b) 达到 最 大 进化 代数 epoch 时 ， 停 止 。 
达到 最 大 进化 代数 之 前 算法 终止 , 则 本 次 运行 结果 为 成 功 。 
实验 种 群 数 目 取 60，epoch 取 50 000。 实 验 均 在 相同 配置 的 计 
算 机 上 ， 使 用 MATLAB 2013b 进行 编程 测试 。 
3.2 测试 结果 分 析 
针对 同一 测试 函数 , 每 种 算法 独立 运行 30 次 , 比较 六 种 算 
法 对 15 个 标准 测试 函数 结果 的 统计 平均 值 (mean )、 方 差 (std)、 
平均 迭代 次 数 (miter〉 及 成 功率 (srate) 。 测 试 结果 如 表 3 所 
示 。 从 表 3 中 可 以 看 出 , 在 30 次 的 独立 运行 中 , 无 论 是 在 单 峰 
测试 函数 、 多 峰 少 极 值 函数 还 是 复杂 的 多 峰 多 极 值 函数 上 ， 
MEBSA 都 能 达到 预 设 精度 ， 且 迭代 次 数 远 远 小 于 其 他 五 种 算 
法 。 充 分 证 明了 MEBSA 不 仅 具 有 较 快 的 收敛 速度 和 强 鲁 棒 性 ， 
还 有 较 强 的 全 局 搜索 性 能 ， 能 够 快速 跳出 局 部 最 优 ， 和 避免 早熟 
收敛 。 


录用 稿 
表 3 六 种 算法 对 15 个 测试 函数 的 测试 结果 
问题 统计 量 ”BSA DS CSO COBSA MABC MEBSA 
mean 3.04e-158 9.12e-257 0 0 5.01e-16 0 
Bl std 9.02e-158 0 0 0 5.83e-17 0 
miter ©50000 50000 23003 10538 50000 2356 
srate 0 0 0.93 1 0 1 
mean 2.22e-95 3.99e-143 0 0 8.51e-16 0 
F2 std 3.70e-95 1.22e-142 0 0 1.10e-16 0 
miter C50000 50000 13851 15496 50000 4045 
srate 0 0 1 1 0 1 
mean 1.56e-158 1.69e-257 0 0 5.19e-16 0 
a std 5.08e-158 0 0 0 6.46e-17 0 
miter ©50000 50000 20375 10452 50000 2158 
srate 0 0 1 1 0 1 
mean 1.44e-96 6.27e-145 0 0 648.6356 0 
F4 std 2.39e-96 2.89e-144 0 0 44.0127 0 
miter ©50000 50000 13845 15347 50000 4082 
srate 0 0 1 1 0 1 
mean 0.0020 0.0086 1.46e+03 0 2.4568et+03 0 
Ee std 0.0014 0.0094 505.9816 0 299.1724 0 
miter ©50000 50000 50000 27097 50000 3924 
srate 0 0 0 1 0 1 
mean 0 0 0 0 8.11e-16 0 
F6 std 0 0 0 0 8.60e-17 0 
miter 48359 38504 6752 10521 50000 2510 
srate 1 1 | 1 0 1 
mean 1.85e-158 5.22e-256 0 0 5.23e-16 0 
二 std 9.33e-158 0 0 0 6.84e-17 0 
miter © $50000 50000 16510 11087 50000 2204 
srate 0 0 0.8 1 0 1 
mean 1.07e-05 2.67e-06 0.85 2.48e-21] 8.87e-16 0 
Fg std 1.12e-05 1.44e-06 2.69 0 1.32e-16 0 
miter ©50000 50000 50000 50000 50000 3825 
srate 0 0 0 0 0 1 
mean 0.6632 0.3699 0.2409 0.0999 2.4199 0 
Fo std 0.1129 0.0535 0.0555 2.3le-l2 0.1375 0 
miter ©50000 50000 50000 50000 50000 2277 
srate 0 0 0 0 0 1 
mean 2.45e-13 9.84e-52 0.0095 0 1.77e-15 0 
BG std 1.32e-12 5.15e-51 0.0153 0 1.47e-15 0 
miter ©50000 50000 30202 15461 50000 4121 
srate 0 0 0.57 1 0 1 
mean 0.0995 0.1658 24.5163 0 0 0 
F11 std 0.3036 0.3771 94.3434 0 0 0 
miter C33252 13523 10757 661 9703 158 
srate 0.5 0.83 0.83 1 1 1 
mean 2.47e-04 0 1.11e-17 0 1.07e-16 0 
i std 0.0014 0 3.39e-17 0 1.07e-16 0 
miter 7379 3864 9803 698 40922 165 
srate 0.97 1 0.90 1 0.33 1 
mean 0 0 0 0 0.9550 0 
E13 std 0 0 0 0 0.0794 0 
miter  é 29357 21685 3039 7938 50000 2930 
srate 1 1 1 ll 0 1 
mean 0 0 2.37e-15 0 0 0 
F14 std 0 0 9.20e-15 0 0 0 
miter 8827 5022 5367 765 3454 170 
srate 1 1 0.93 1 1 1 
mean 5.30e-14 0 1.89e-14 0 3.79e-15 0 
i std 6.49e-14 0 5.24e-14 0 2.07e-14 0 
miter 31886 6438 9656 652 15839 162 
srate 0.57 1 0.87 1 0.97 1 
为 了 进一步 直观 说 明 MEBSA 算法 的 先进 性 ， 图 1~4 给 出 
了 典型 的 四 种 测试 函数 的 收敛 曲线 。 其 中 横 坐 标 为 迭代 次 数 ， 
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纵 坐 标 为 函数 值 。( 注 :为 了 方便 进化 上 


函数 值 取 以 10 为 底 旧 


C "INaxX IV Wa 
等 : 基于 变异 交叉 方程 与 进化 选择 机 制 的 回溯 优化 改进 算 ; 


人 + 口 工 | 
Pl -|- 
| 


H 线 的 显示 和 观察 ,本 文 对 


函数 值 


函数 值 


函数 值 


函数 值 


的 对 数 ) 。 
F3 
-BSA 
-DS 
| -CSO 
-COBSA 
-MABC 
| 一 一 MEBSA 
2 3 4 5 
代数 4 
x10 
图 1 F3 函数 的 收敛 曲线 
F6 
-BSA 
-DS 
| -CSO 
-COBSA 
-MABC 
. -一 一 MEBSA 
2 3 4 5 
人 x 10 
图 2 F6 函数 的 收敛 曲线 
F8 
-BSA 
~ :DS 
~ -CSO 
se 一 'COBSA 
Se -MABC 
~~、. 一 一 MEBSA 
2 3 4 5 
代数 本 104 
图 3 F8 函数 的 收敛 曲线 
F10 
本 a -BSA 
-DS 
-CSO 
-COBSA 
-MABC 
| -一 一 MEBSA 
和 
AN 
由 
2 3 4 5 
代数 
x 10 
图 4 Fl10 函数 的 收敛 曲线 
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表 可 以 看 出 ， 针 对 F3、F6、F10 等 其 他 五 种 算法 也 
能 部 分 成 功 运行 的 函数 ，MEBSA 能 以 最 少 的 迭代 次 数 以 及 最 
快 的 收敛 速度 达到 预 设 精度 ; 针对 F8、F9 等 其 他 五 种 算法 均 
不 能 成 功 的 函数 ，MEBSA 依然 能 以 较 快 速度 成 功 收敛 。 进 一 
步 展 示 了 MEBSA 算法 的 优越 性 


o 
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为 了 解决 标准 BSA 算法 收敛 速度 相对 较 慢 及 搜索 精度 不 
高 的 问题 ， 本 文 提出 了 新 的 变异 尺度 系数 F， 并 加 入 最 优 个 体 
的 引导 ， 构 造 了 一 个 新 的 交叉 方程 ， 同 时 ， 为 了 防止 算法 陷入 
局 部 最 优 ， 提 出 了 进化 选择 机 制 ， 引 入 差分 进化 变异 算 子 。 最 
后 ,选取 15 个 测试 函数 , 并 与 五 种 算法 进行 了 对 比 。 结 果 表明 ， 
本 文 改 进 算法 MEBSA 有 效 提高 了 算法 的 各 项 搜索 性 能 。 


参考 文献 : 


[1] Holland J H. Adaptation in natural and artificial systems [M]. [S. 1.] : Ann 
Arbor. University of Michigan Press, 1975. 

[2] Eberhart R, Kennedy J. A new optimizer using particle Swarm theory [C1// 
Proc of Micro Machine and Human Science. Washington DC: IEEE 
Computer Society, 1995: 39-43. 

[3] Storn R, Price K. Differential evolution: a simple and efficient heuristic for 
global optimization over continuous spaces [J]. Journal of Global 
Optimization, 1997, 11 (4): 341-359. 

[4] Karboga D, Basturk B. A powerful and efficient algorithm for numerical 
function optimization: artificial bee colony (ABC) algorithm [J]. Journal of 
Global Optimization, 2007, 39 (3): 459-471. 

[5] Zhu Guopu, Kwong S. Gbest-guided artificial bee colony algorithm for 
numerical function optimization [J]. Applied Mathematics & Computation, 


2010, 217 (7): 3166-3173. 


Chi naXiv 合 作 期刊 
赵 琳 敬 ， 等 : 芝 下 交 关 和 择机 制 的 杏 湖 优化 改进 法 


[6] 李 文 斌 , 陈 变 瑛 , 贺 裔 朝 , 等 . 邻 域 结 构 为 复杂 网 络 的 差分 演化 算法 
四 . 计算 机 应 用 研究 , 2016, 33 (2): 370-374. 

[7] 姜 凤 利 ,， 张 宇 ， 王 永 刚 ， 等 . 一 种 基于 引导 策略 的 自 适 应 粒子 群 优化 算 
法 [J]. 计算 机 应 用 研究 ,2017, 34 (12): 3599-3602 

[8] Meng Xianbing, Liu Yu, Gao Xiaozhi, et al. A new bio-inspired algorithm: 
chicken Swarm optimization [C]/ Proc of International Conference on 
Swarm Intelligence. Hefei: Springer International Publishing, 2014: 86-94. 

[9] Zheng Yujun. Water wave optimization: a new nature-inspired metaheuristic 
[J]. Computers and Operations Research, 2015, 55: 1-11. 

[10] 高 江 锦 , 杨 烙 . 免疫 进化 否定 选择 算法 [J]. 计算 机 应 用 研究 , 2017, 34 
(5): 1293-1297. 


[11 


Civicioglu P. Backtracking search optimization algorithm for numerical 
optimization problems [J]. Applied Mathematics and Computation, 2013, 
219 (15): 8121-8144. 

[12] 田 文 饥 , 刘 三 阳 , 王 晓 娟 . 基于 差分 进化 的 回溯 搜索 优化 算法 研究 与 改 
进 [J]. 计算 机 应 用 研究 , 2015, 32 (6): 1653-1656. 

[13] Zhao Lei, Jia Zhicheng, Chen Lei, et al. Improved backtracking search 
algorithm based on population control factor and optimal learning strategy 
[J]. Mathematical Problems in Engineering, 2017, 5 (1): 26-36. 

[14] Civicioglu P. Transforming geocentric cartesian coordinates to geodetic 
coordinates by using differential search algorithm [J]. Computers and 
Geosciences, 2012, 46 (3): 229-247. 

[15] Gao Weifeng, Liu Sanyang. A modified artificial bee colony algorithm [J]. 
Computers and Operations Research, 2012, 39 (3): 687-697. 

[16] Liang Jing, Qin A K, Suganthan P N, et al. Comprehensive learning particle 
Swarm optimizer for global optimization of multimodal functions [J]. IEEE 
Trans on Evolutionary Computation, 2006, 10 (3): 281-295. 

[17] Karaboga D, Akay B. A comparative study of artificial bee colony algorithm 


[J]. Applied Mathematics and Computation, 2009, 214 (1): 108-132. 


